Classical and quantum anomalous diffusion in a system of 25-kicked Quantum Rotors 
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We study the dynamics of cold atoms subjected to pairs of closely time-spaced 5-kicks from 
standing waves of light. The classical phase space of this system is partitioned into momentum cells 
separated by trapping regions. In a certain range of parameters it is shown that the classical motion 
is well described by a process of anomalous diffusion. We investigate in detail the impact of the 
underlying classical anomalous diffusion on the quantum dynamics with special emphasis on the 
phenomenon of dynamical localization. Based on the study of the quantum density of probability, 
its second moment and the return probability we identify a region of weak dynamical localization 
where the quantum diffusion is still anomalous but the diffusion rate is slower than in the classical 
case. Moreover we examine how other relevant time scales such as the quantum-classical breaking 
time or the one related to the beginning of full dynamical localization are modified by the classical 
anomalous diffusion. Finally we discuss the relevance of our results for the understanding of the 
role of classical cantori in quantum mechanics. 

PACS numbers: 32.80.Pj, 05.45.Mt, 05.60.-k 



I. INTRODUCTION 

The effect of the underlying classical dynamics on the 
quantum motion has been a recurrent topic of research 
since the early days of the quantum theory. In recent 
years the use of experimental techniques based on ultra- 
cold atoms in optical lattices [l| has permitted to study 
in great detail the role of classical mechanics in simple 
quantum systems. In these experiments a very dilute al- 
most free gas of atoms (Cs and Rb) is cooled down to 
temperatures of the order of tens and then interacts 
with an optical lattice. In its simplest form, the opti- 
cal lattice consists of two laser beams prepared in such a 
way that the resulting interference pattern is a stationary 
plane wave in space. The laser frequency is tuned close 
to a resonance of the atomic system in order to enhance 
the atom-laser coupling but not too close to avoid spon- 
taneous emission. In this limit the laser-atom system 
can be considered as a point particle in a sine potential, 
namely, the quantum pendulum. If the laser is turned 
on and off in a series of short periodic pulses, the result- 
ing system is very well approximated by the so called 
quantum kicked rotor (QKR) (see Ref. W\ for a review) 
extensively studied in the context of quantum chaos. 
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Kcos{q)^5{t-Tn). 



(1.1) 



For short time scales, quantum and classical motion 
agrees. However quantum diffusion is eventually sup- 
pressed due to destructive interference that localize 
eigenstates in momentum space. This counter-intuitive 
feature, usually referred to as dynamical localization 3], 
was fully understood [3| after mapping the kicked rotor 
problem onto a short range one dimensional disordered 



system where localization is well established. The theo- 
retical predictions of Ref. [3| were eventually confirmed 
experimentally [l| (see also Ref. j3]) by using the cold 
atoms techniques mentioned previously. The standard 
kicked rotor is thus an ideal candidate for the study of the 
quantum properties of one dimensional systems whose 
classical motion is diffusive. A natural question to ask is 
whether this analysis can be extended to other types of 
(anomalous) diffusive motion. 

For values of the kick strength K in Eq. (|l.ip suffi- 
ciently small, the classical phase space is composed of 
chaotic and integrable parts and, for certain initial condi- 
tions, the classical motion is well described by a process 
of anomalous diffusion. The quantum transport prop- 
erties in systems with a mixed phase space [6(] depend 
strongly on the details of the Hamiltonian ?|. Even for 
a given configuration many types of anomalous diffusion 
are observed depending on the initial conditions or the 
time scales studied Q. This lack of universality makes 
it difhcult to precisely assess the effect of the underlying 
classical anomalous diffusion on the quantum dynamics. 

The situation is different if the smooth sinusoidal op- 
tical potential is replaced by a potential with a loga- 
rithmic or power-law singularity [9| . The classical phase 
space is homogeneous but still the classical motion is su- 
perdiffusive. In the quantum realm, as a consequence 
of interference effects, the particle still diffuses but at a 
slower rate. In fact, for certain types of singularities, full 
dynamical localization never occurs and diffusion per- 
sists at all times. In other cases exponential localization 
is eventually observed but anomalous diffusion different 
from the classical one is still observed for shorter times. 
The classical density of probability P{p, t) in this region 
is accurately described by the solution of an anomalous 
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Fokker-Planck equation [T3|- We note that, unlike the 
case of normal diffusion, the information obtained from 
the knowledge of a few moments of the density of proba- 
bility may not be sufficient to fully characterize the clas- 
sical motion (lo| . Thus for a correct understanding of 
these systems it is essential a detailed investigation of 

The models studied in Ref. @ are non-KAM but clas- 
sical anomalous diffusion can also exist in KAM systems. 
One example is the kicked rotor with a smooth potential 
but subjected to pairs of cl osely time-spaced kicks: the 
2(5-kicked rotor (2(5-KR) [H . 

The dynamics of this model in a certain region of pa- 
rameters has already been investigated in the literature: 
theoretically in Ref. [ll| and experimentally in Ref. [l^] . 
Unlike the single kicked rotor the classical dynamics is 
strongly correlated. The momentum space is divided into 
regions of fast momentum diffusion separated by porous 
boundaries, ie narrow trapping regions where classical 
trajectories 'stick' for relatively long periods. The 26- 
KR trajectory spends considerable time trapped in a cell 
before escaping to the next. 

In this paper we provide a detailed account of the type 
of classical and quantum anomalous diffusion associated 
with this motion. We shall restrict ourselves to the one 
cell region, namely, to time/momentum scales such that 
a particle initially trapped manages to escape and even- 
tually reaches a new trapping region. We aim a bet- 
ter understanding of how generic features of the classical 
anomalous diffusion affect the quantum motion in a KAM 
system. We restrict ourselves to a region of parameters 
such that the dynamics is generic, namely, the classical 
phase space is fully chaotic with no island of stability. 
However the trapping-leaking mechanism of our model 
still causes strong deviations from the single kicked rotor 
results. Other important motivation to study the mo- 
tion of a fully chaotic 26-KR is that it mimics the effect 
of cantori in chaotic systems. Thus a particle typically 
gets trapped in a cantorus for a long time until it escapes 
to the chaotic sea. We argue the findings of this paper 
shed light about the role of classical cantori in quantum 
mechanics [isl. [T5l|. 

The organization of the paper is as follows: in the next 
section we introduce the model, review some of its more 
relevant dynamical features and discuss the region of pa- 
rameters to be studied. In section III, we study the classi- 
cal and quantum transport properties. Among others we 
analyze the classical and quantum density of probability, 
the classical-quantum breaking time and the return prob- 
ability. We aim to describe how dynamical localization 
arises in systems whose classical diffusion is anomalous 
and how other relevant scales of the problem such as the 
quantum-classical breaking time are affected by the un- 
derlying classical dynamics. 

In summary, our main new results are: 1) In the re- 
gion of interest both classical and quantum dynamics is 
well characterized by a process of anomalous diffusion. 
2) We identify two routes to dynamical localization as 



a function of H: for h > he {he is function of the kick 
strength and the separation between pairs of kicks) stan- 
dard dynamical localization occurs in the trapping region 
and the particle does not escape; for h < he the cen- 
tral part probability density is exponential and almost 
time independent. However diffusion does not stop since 
eventually the particle escapes from the trapping region. 
True dynamical localization in this case occurs for time 
scales much longer than the typical time to escape from 
the trapping region. For intermediate times we find the 
quantum diffusion is anomalous but slower than the clas- 
sical one. 3) The classical anomalous diffusion induces a 
fractional scaling with h in different quantities of interest 
such as the quantum-classical breaking time. 4) The 2(5- 
KR can be used as a simplified model to study the effect 
of cantori in classical and quantum mechanics. 

II. THE MODEL 

We consider a system with a Hamiltonian correspond- 
ing to a sequence of closely spaced pairs of kicks: 

2 

^-Kcosx^[S{t-nT) + S{t-nT + e)], 

n 

where e ^ T is a short time interval and K is the kick- 
strength. 

A. Classical dynamics 

The classical map for the 26-KP is a straightforward 
extension of the Standard Map: 

Pn+l = P7i + KsinX„; Pn+2 = Pn+i + Ksinxn+i 
Xn+l ^ Xn+ ePn+i; Xn+2 = Xn+l + TPn+2 

(2.1) 

where e is a very short time interval between two kicks 
in a pair and r = T — e is the (much longer) time in- 
terval between the pairs. Clearly, the limit e = r or 
corresponds to the Standard Map, which describes the 
classical dynamics of the quantum kicked rotor: 

Pi+i = pt + Ksinxt, 

x.i+1 = Xi +Pi+i. (2.2) 

Particles with momenta pq ~ (2m -I- l)n/e and m = 
0, ±1, ±2, • • • (relative to the optical lattice) are confined 
in momentum trapping regions and absorb little energy; 
conversely, particles prepared near po ~ 2rmr/e experi- 
ence rapid energy growth up to localization. The ba- 
sic mechanism of trapping is fairly intuitive: atoms for 
which Po — {2m + l)7r/e experience an impulse sina; 
followed by another one ~ Ksm{x + n) which in effect 
cancels the first. Over time, however, there is a gradual 
de-phasing of this classical 'anti- resonant' process. A the- 
oretical study of the classical diffusion in this system 
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found anomalous momentum diffusion for any pq and in- 
termediate time scales. It was also observed long-ranged 
corrections to the uncorrelated diffusion rate not present 
in the standard kicked rotor. Inside the trapping region 
diffusion is normal but the coefficient of diffusion D has 
a dependence on K a,s D ^ , similar to the one found 
in single kicked rotors in a region of the classical phase 
space densely populated by cantori .13 .]. By contrast, if 
the phase space is fully chaotic it is expected D ^ K^. 
This suggests that our model may reproduce to a good 
approximation the effect of cantori in a KAM system. 

The size in momentum space of the trapping region 
Sp strongly depends on the parameters e, K defining the 
model. 

For an accurate analysis of the trapped region is neces- 
sary that: i) the particle is initially trapped, ii) the par- 
ticle dwells on average a time long enough in this region 
before it escapes from it. In Ref. it was concluded 
that this implies the criterion Ke <C 1. However, if Ke is 
too small, the phase space would be too regular. Since we 
are interested in generic feature of the motion, we take 
iii) T ^ 1. This implies less correlations between suc- 
cessive space positions of the particle. We thus restrict 
ourselves through the paper to the range Ke <C 1 and 
T 3> 1 with initial conditions inside the trapping region. 

B. Quantum dynamics 

The time evolution operator for this system can be 
written as 

= e^'^e*"^ ''"''"'e"*^ e^f (2.3) 
In a basis of plane waves, Lf^ has matrix elements 

TT'^ — Tj}^<^'^ Tj2-kick _ ^-i^^^ -l-m 

Y,Ji~kiKn) Jk-m{Kn) e-'^ (2.4) 
fc 

where Kn — K/h and Jn{Kh) is integer Bessel functions 
of the first kind. It is easy to see that Uf^^'^'^^ is in- 
variant if the products K^^ = Ke and = he are kept 

constant; while the free propagator ul^'^'' — e^*'^~^^ sim- 
ply contributes a near-random phase. Thus the results 
are quite insensitive to the magnitude of r = T — e pro- 
vided that T » 1. We will stick to Ke <^ 1 and r 1 in 
all numerical calculations. 

The result in Eq. (j2.4|) may be compared with the 
one-kick map in Eq. (|2.ip 

f^/™ =6-^"^ Jl-ra{Kn). (2.5) 

The one-kick matrix for the QKR has a well-studied 
band-structure: since Ji_„i{x) ~ for \l — m\ ^ x, we 
can define a bandwidth for C/(°); i.e. h = Kf^ (this is 
strictly a /ia//-bandwidth) which is independent of the 



angular momenta I and m. However, this is not the case 
for the matrix of V^. 

Assuming \l — m\ is small it was shown in Ref. [Tlj 
that 

t/L « e-''* Ji-™ [2Kh cos {lh,/2)] (2.6) 

where the phase $ = (PT -I- elm + eP)h/2 + 7r(/ - 
to)/2. Hence we infer a momentum dependent band- 
width, b{p) = 2KfiCOs{pe/2). 

While [/'■*'•' has a constant bandwidth, the bandwidth 
for the matrix of [/' oscillates with I from a maximum 
value bmax = 2Kfi, equivalent to twice the bandwidth of 
U^^\ down to a minimum value 6mm ~ 0. is thus 
partitioned into sub-matrices of dimension N = 27: /h^ 
corresponding precisely to the momentum cells of width 
Ap — Nh observed in experiments . 

In this paper we aim to study both the evolution of an 
initially given wave-packet and its low order moments. 
The structure of the evolution operator [/^ allows an ef- 
ficient and accurate numerical calculation of these quan- 
tities. The action of the operator If^ on a quantum state 
can be decomposed into four steps: two associated with 
the kicks, which are diagonal in the space presentation; 
and two associated with the free rotations between neigh- 
boring kicks, which are diagonal in the momentum pre- 
sentation. We thus take the space and momentum rep- 
resentations alternatively to facilitate the calculations. 
The transformation between the representations is effi- 
ciently carried out with the fast Fourier transformation 
(FFT) algorithm. In our numerical simulations we can 
set the size of the basis up to 2^° which is big enough to 
guarantee a double precision accuracy. 



III. RESULTS: ANOMALOUS DIFFUSION AND 
DYNAMICAL LOCALIZATION 

In this section we investigate the classical and quan- 
tum dynamics of the 2(5-KR. Our main motivation is to 
describe the effect on the quantum motion of generic 
dynamical features of the classical dynamics such as 
the trapping-escaping mechanism. We mainly focus on 
time/momentum scales such that the motion is con- 
fined inside one cell. Initial conditions are always chosen 
within the trapping region. For a study in the multi-cell 
region we refer to Ref. [HI, . 

As was mentioned previously, we restrict ourselves to 
the window of parameters Ke ^ 1. In addition it is also 
imposed r ^ 1 in order to remove correlations that make 
the dynamics non generic. In the numerical calculations 
T, K, and e are fixed within the above limits. Our main 
conclusions do not depend on the specific value of the 
parameters. Quantum effects are investigated by varying 

h. 
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FIG. 1: (Color online) Comparison of quantum and classi- 
cal energy diffusion (from bottom to top: ft = 4 x 10~^, 1 x 
10~^, 5 X 10~*, 2.5 X 10~* and classical respectively) for K = 
0.2, e — 0.25 and r = 10*. The two figures only differ in the 
scale of time/momentum represented. Fig. (a) shows the effect 
of leaking from the trapping region in longer time scales while 
Fig.(b) describes diffusion inside the trapping region. The 
inset of Fig.(b) shows the dependence of quantum-classical 
breaking time tc on h; the best fitting (thin solid line) sug- 
gests that tc ~ ;i-i-39±o.o5 ^ ^ -pj^j^ contrast 
with the single kicked rotor where tc ~ 



A. Energy diffusion 

The classical dynamics in the trapping region was in- 
vestigated analytically in Ref . [13 • It was found that en- 
ergy diffusion (p^(i)) ~ Dt increases linearly with time. 
However, the dependence of D ~ on the kicking 
strength is different from the single kicked rotor predic- 
tion D in the chaotic region but similar to the one 
observed if the classical phase space is populated by can- 
tori This suggests that our model captures generic 
features of the effect of cantori in classical mechanics. 

Eventually the particle reaches the edge of the trap- 
ping region and leaks outside. Diffusion becomes then 
anomalous until the particle gets trapped again. As is 
shown in Fig. 1, the parameters K^e and r were cho- 
sen such that the typical time for considerable leaking 
is tieak ~ 200 3> 1. Before this time the leakage is still 



possible but with a negligible probability. As the parti- 
cle approaches the outer boundaries of the two cells that 
sandwich the trapping region (at t « 6 x lO"^), the dif- 
fusion begin to slow down. The time unit we take in all 
figures is the number of the kick pairs. 

In the quantum case the dynamics is richer. If we re- 
strict ourselves to the one cell region three time scales are 
distinguished. In a first stage, t < t^, the quantum av- 
eraged energy agrees with its classical counterpart. As a 
result of the peculiarities of the classical dynamics, scal- 
ing with h of the classical-quantum breaking time tc may 
be different from tc ~ the result for a single kicked 
rotor. For longer times, tc < t < td, interference effects 
start to be important. As a consequence we expect the 
particle still diffuses, {p'^{t)) « Dquant^ , but at a rate 
which decreases as h increases and it is lower than the 
classical one; namely 7 < 1 and Dquan < Ddas- This is 
a novel regime caused by the interplay of destructive in- 
terference and classical anomalous diffusion. This stage 
lasts up to a time id in which full dynamical localization 
occurs and diffusion is totally arrested. 

Our numerical results confirm the above qualitative 
picture. Results for quantum and classical energy diffu- 
sion (p^lt)) (p = hk for the quantum case) as a function 
of time t for T = 10^ and Ke = 0.05 < 1 for different ft's 
are shown in Fig. 1. In the classical case 2 x 10^ initial 
conditions uniformly distributed along = 7r/e (located 
at the center of the trapping region) have been utilized 
for the ensemble average. The quantum initial condition 
is set as the momentum eigenstate |^o) with the integer 
closest to po I h. 

We restrict ourselves to time/momentum scales such 
that the motion is confined to one cell. We observe that 
for short times both classical and quantum results coin- 
cide up to a certain breaking time tc which has a frac- 
tional scaling with h, different from the single kick rotor 
case (see insert of Fig. lb, where tc is evaluated as the 
maximum time such that the deviation between classi- 
cal and quantum {p'^{t)) is below 20%). For later times 
the quantum particle still diffuses normally but with a 
smaller diffusion coefficient Dquan (see Tab. I). This is 
a direct consequence of destructive interference effects. 
We relate this new region of weak dynamical localiza- 
tion to the effect of classical trapping on the quantum 
dynamics. This interpretation is reinforced by the ob- 
served dependence of Dquan on h in the trapping region 
(see Tab. I). For h sufficiently small [h < 10"*), Dquan is 
close to the classical prediction. For larger h we observe 
a growing dependence on h though the exponent 7 ~ 1 
is not modified. It decreases as h increases. This situa- 
tion lasts until a maximum h, denoted as Hmax , such that 
the breaking-time tc coincides with the time td in which 
full dynamical localization starts. For the parameters 
utilized hmax ~ 6 X 10^'^. 

We have observed an additional h scale, denoted as he, 
relevant for a precise characterization of the trapping- 
leaking mechanism. It is defined by the smallest h such 
that td < ticak- For the parameters used he ~ < 
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TABLE L Breaking time tc, dynamical localization time td, 
and classical/quantum diffusion coefficient against h for K — 
0.2, e = 0.25, and r = 10*. D^ias « 1.83 x 10~^ which is 
evaluated over t < tieak ~ 200. Dquan is evaluated over the 
time range tc < t < tieak for h < he ~ 10~^ and tc < t < td 
otherwise. 



Sma:!;- Foi He < H < Hmax, td < iieafe, -D,„a« IS cloSC tO 

zero, and the wave-packet is well localized in the trapping 
region. In this range as h approaches hmax, td tends to 
tc- For h < he, the time scale td related to the dynami- 
cal localization increases dramatically. However quantum 
energy diffusion is still linear in time (like the classical 
one) up to tieak- After t « t;eafc the classical particle 
starts to leave the trapping region and the classical mo- 
tion become superdiffusive. The quantum dynamics for 
h < he is also superdiffusive, but with a smaller leaking 
rate. The smaller the h, the closer the leaking rate to the 
classical one. Eventually w 6 x 10"^) a new trapping 
region is approached. As a consequence, both classical 
and quantum motions are slowed down again. 

We note that the region of weak dynamical localiza- 
tion characterized by a linear time dependence of the en- 
ergy evolution is not present in the standard kicked rotor 
where there is no intermediate region between classical 
diffusion and full quantum dynamical localization. This 
is an indication that even though the quantum diffusion 
is still normal inside the trapping region the 2(5-KR is 
essentially different from the single kicked rotor. 

For longer times t > td diffusion stops as a consequence 
of standard dynamical localization similar to the one ob- 
served in the single quantum kicked rotor. Our results 
thus suggest that in order to observe genuine quantum 
anomalous diffusion the value of h must be such that 
td ^ tc- This corresponds with h < he- We remark that 
this condition should be met by any experiment aiming 
to confirm the results reported in this paper. 

Although not shown, we have confirmed that the de- 
pendence of the classical diffusion constant is £> ~ 
similar to the case of a single kick rotor in a region popu- 
lated by classical cantori. This together with the anoma- 
lous dependence with h also found in studies of the role of 
cantori in quantum mechanics is a further indication 
that the 2(5— KR in the region of parameters studied in 
this paper can be utilized as a effective model to investi- 
gate the role of classical cantori in quantum mechanics. 



Finally we note that ^ tis only a necessary con- 

dition for normal diffusion but this by no means assures 
that the density of probability is Gaussian- like [l3| • In- 
deed, in the next section we will show that the quantum 
and classical density of probability of our model strongly 
deviate from the normal diffusion prediction. 

B. Density of probability 

The density of probability of finding a particle with 
momentum p after a time t from a given initial state 

= \lo) is given by P,(p,i) = Pq{k,t) = 
with p = kh- The parameter set K = 0.2, e = 0.05 and 
T — 10* chosen permits us to study generic features of 
the motion. By generic we mean the trapping region is 
large enough to be studied and the classical phase space 
has not stability island; namely, Ke ^1. In addition 
r ^ 1, so consecutive pairs of kicks are uncorrelated. In 
all cases our initial state is located in the trapping region. 

1. Classical density of probability 

The classical P{p, t) is obtained by evolving the classi- 
cal equation of motion for 2 x 10^ different initial condi- 
tions at the center of the trapping region. Positions are 
uniformly distributed along the interval (— 7r,7r). 

P{p, t) is approximated as the set of the states of the 
whole ensemble that falls in (p— Ap/2,p+Ap/2) at time t. 
In the numerical calculations we set Ap = 0.02V(p2(i)}. 

We distinguish the following regions in the classical 
density of probability (see Fig. 2a): For short times such 
that the particle is well inside the trapping region the 
diffusion is normal and P{p, t) is Gaussian (see inset of 
Fig. 2a). As the boundary of the trapping region is 
approached we observe a gradual crossover from normal 
to anomalous (super) diffusion. For small momentum 
P{p, t) is still Gaussian as leaking to the outside region 
is weak. As time approaches ti^ak ~ 200, the typical 
time to reach the edge of the trapping region, the cen- 
tral (small momentum) Gaussian region becomes smaller 
and smaller. Meanwhile, the outskirts bend down and 
a power-like behavior typical of anomalous diffusion is 
observed. 

For longer times, t > tieak, P(p, t) ~ p^" with a a de- 
creasing function of time (Fig. 2a). For t oo, a 0. 
This behavior is not surprising due to the transient na- 
ture of the trapping-leaking mechanism. The power-law 
decay is a direct consequence of the enhanced diffusion 
once the particle has escaped from the trapped region. 
However, precisely due to this fast diffusion, the particle 
soon approaches a new trapped region where diffusion is 
slowed down again until the particle leaks from the trap. 
As a consequence of this bottle neck, the exponent a de- 
creases and the density of probability in the region of 
fast diffusion gradually increases with time. In addition 
a sharp jump in P{p, t) is observed between the enhanced 
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diffusion zone and the new trapping region (see Fig. 2a 
for t = 6000 where a steep drop of P{p, t) happens at 
around |p — po| — 27r/e corresponding to the new trap- 
ping region). For long times t — > oo, P(p,t) resembles a 
staircase, flat between trapping regions and discontinu- 
ous in the trapping areas. This is a simple consequence 
of the fact that the dwelling time in the trapping region 
is typically much longer than the one needed to travel 
between two consecutive traps. We shall see in the next 
section that quantum mechanically the situation is dif- 
ferent. 



2. Quantum density of probability 
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FIG. 2: (Color online) Distribution of the classical (a) and 
quantum probability density with ft « 5 x 10~* (b) and 
3 X 10~^ (c) respectively. Classical P{p, t) is characterized 
by a Gaussian profile (inset of Fig. (a)) before leaking 
from the trap occurs {t < ti^ak ~ 200). For later times 
tieak <t <td, P{p,t) devclops powcr-law tails, P{p,t) ~ p~" , 
typical of anomalous diffusion. The exponent a varies in time 
{a 4.5, 2.9, 2.3, 1.3 for t = 500, 1000, 2000 and 6000 respec- 
tively). In the quantum case for small ft leaking leads to 
power- law tails as well (see Fig. (b)). The best fitting ex- 
ponent, indicated by thin black lines, is a ~ 7.1,4.2,3.0,0.94 
for t = 500,1000,2000 and 6000 respectively). For larger ft 
(see Fig. (c)), dynamical localization suppresses the diffusion 
before leaking occurs. K = 0.2, e = 0.25 and r = 10*. 



The quantum density of probability Pq {p, t) was cal- 
culated from an initial condition |'0(O)) — \lo) with 
^0 = T^/fie- For given values of h and e the right hand 
side may not be an integer. Then the nearest value of 
h was chosen to ensure this condition. Thus po = 
exactly for both classical and quantum calculations. As 
in the classical case, Pq{p,t) was evaluated by summing 
up the probability of falling in a bin of width Ap. We set 
Ap = 0.0315 for /i « 5 X lO"" (Fig. 2b) and Ap = 0.006 
for h ^ 3 X 10~^ (Fig. 2c). Hence it is in fact a coarse- 
grained result where part of quantum fluctuations have 
been suppressed. It is more instructive to start our anal- 
ysis of the quantum probability density with a general 
account of the expected effect of the classical motion on 
the quantum dynamics. As in the single kicked rotor we 
expect to observe dynamical localization for sufhciently 
long times t > td However the trapping (and even- 
tual release) mechanism changes qualitatively the route 
to dynamical localization in our model. In fact we distin- 
guish two different scales for full dynamical localization. 
For h > he, full dynamical localization will occur well 
inside the trapping region where the classical diffusion is 
normal. As a consequence the particle will not have time 
to escape, {td < i/eafe ~ 200 with our parameters) and the 
density of probability will be time independent and will 
decay exponentially as a function of the momentum p. 
These are typical signatures of exponential localization 
similar to the one observed in a single kicked rotor. The 
situation is different if h < he is small enough such that 
leaking to the enhanced diffusion region is possible. In 
this case we may still observe typical features of dynam- 
ical localization within the trapping region. However for 
larger momenta the density of probability develops time- 
dependent power law tails typical of anomalous diffusion. 
Full dynamical localization will not typically take place 
until the next trapping region is reached or later if h is 
small enough. This revival of quantum diffusion is a typ- 
ical feature of our model that it is not observed in the 
single kicked rotor. Below we provide a more detailed ac- 
count of the relevant time and momentum scales (see Fig. 
2b) for an accurate description of the quantum density 
of probability: 

1. For t,p well within the trapping region: t < ti^ak 
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and |p— Pol < Pieak ~ 0.1. After a narrow region in which 
classical and quantum results fully agree (see inset Fig. 
2a), we observe in the quantum case a cusp for p ^ po 
caused by the incipient dynamical localization in the core 
of the trapping region. However for larger momentum 
the distribution is Gaussian in agreement with previous 
results (see Fig. 1) for the energy diffusion. 

2. For p within the trapping region (|p — po| < Pieak) 
but t > tieak- The core of the quantum probability is 
still exponentially localized (see Fig. 2c) but the out- 
skirts Pq{p, t) ~ develops a power-law tail (see Fig. 
2b) typical of anomalous diffusion. In this region the 
exponent a depends on h but it is time independent. 

3. For p outside or close to the edge of the trap- 
ping region and t > tieak- The quantum probability 
Pq {p, t) 1 /p" develops power-law tails but in this case 
the exponent a decreases with time and eventually tends 
to zero for t — > cx) (see Fig. 2b). This is in principle 
surprising. It is in general expected that destructive in- 
terference slows down the quantum motion. As a conse- 
quence the quantum density of probability should decay 
faster than the classical one. However there is a simple 
explanation for this behavior: in a first stage, the quan- 
tum exponent a is larger than the classical one as a conse- 
quence of destructive interference effects. For later times, 
the classical probability between two trapping region be- 
come smaller than the quantum one since the classical 
particle leaks faster into the next cell. At the same time 
destructive interference slows down the quantum motion 
in this region. Eventually Pq{p,t) saturates, a — > be- 
fore the next trapping region is reached (see t = 6000 in 
Fig. 2a and Fig. 2b). 

4. For times sufficiently long t > quantum destruc- 
tive interference effects dominate, standard dynamical 
localization takes place, and diffusion stops. Pq(p,t) be- 
comes time independent and decays exponentially as a 
function of the momentum. As was mentioned previ- 
ously the value of td depends dramatically on whether 
dynamical localization occurs in the trapping region or 
after it has escaped from the trapping region. In the for- 
mer case td < tieak ~ 200 and in the latter td > 10^ with 
no intermediate values. 



C. Return Probability 

Having studied diffusion in momentum for a fixed time 
t we now look at the explicit dependence in time. In 
order to proceed we calculate return probabilities P{t) 
of a wave-packet as a function of time. We average over 
N initial starting conditions close to the center of the 
trapping region to suppress quantum fluctuations, 

P{t)^il/N)f2{Piit)), (3.1) 
h 

where Pi(t) = |(V'(t)|V'!(i = 0))^. The initial condition 
was taken to be an angular momentum eigenstate \^i(t = 




FIG. 3: (Color online) Return probability as a function of 
time for K — 5, e = 0.04 and r = 10"*. We present results 
for different values of h together with the classical prediction. 
In the trapping region {{t < tieak ~ 30) in this case) the 
return probability decays as a power-law P(t) ~ t The 
exponent 7, given by the best linear fitting (indicated by a 
thin black line), decreases as h increases. We have obtained 
7 « 0.50, 0.46, 0.40 and 0.28 for the classical and quantum 
cases with h = 2~^,2~'' and 2~^ respectively. 



0)) = 1^). The results were further averaged from li « 
n/h^ — N/2 to h = h + N. In order to ensure the initial 
conditions are within the trapping region, the value N is 
set to be iV « Sp/h with Sp the width of the trapping 
region. 

The return probability provides valuable information 
about the degree of localization of a system. Indeed it 
was already utilized in the landmark paper by Anderson 
fv^ about localization. A non zero P{t) in the t 00 
limit is a signature of exponential localization. On the 
other hand P{t) oc t"''/^ with d the spatial dimensional- 
ity is typical of fully delocalized eigenstates (normal dif- 
fusion). In the theory of localization a power-law decay 
P{t) cx l/t"' with 7 < d/2 is a signature of localization 
typical of a disordered conductor close to a localization 
transition. Such slow decay has also been related to the 
effect of cantori in transport [18] properties. 

Our results are summarized in Fig. 3. For the sake 
of comparison, we first study the classical counterpart of 
P{t). Again 2 x 10*" initial states uniformly located in 
space and with po = ir/e were evolved. P(t) was eval- 
uated as the set of states that fall in the momentum 
region |p — po| < 0.004 at time t. For the parameters 
{K = 5, e = 0.04 and r = 10*) chosen for the calculation, 
the width of the trapping region Sp w 4, and tieak ~ 30. 
For t < tieak, P{t) ^ t^'^'^, in agreement with the Gaus- 
sian probability density distribution obtained in the pre- 
vious section. For t > tieak , leaking gradually takes over 
and after a crossover P{t) undergoes a faster power-law 
decay, P{t) ~ t~^-^ until the next trap is reached. 

Within the trapping region the quantum P{t) decays 
as t~'^ with 7 < 0.5, a decreasing function of h. The fact 
that 7 < d/2 is a clear indication that destructive inter- 
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ference is already slowing down the quantum diffusion. 
Not surprisingly, 7 approaches its classical limit, 0.5, as 

For sufficiently large ft > 0.1 the particle cannot es- 
capes from the trapping region and P{t) becomes con- 
stant (see for example P{t) for h = 2^^ in Fig. 3). This 
is fully consistent with previous results for the energy 
diffusion and the density of probability. 

For smaller h the particle escapes from the trapping re- 
gion and full dynamical localization occurs only for times 
much longer than those shown in Fig. 3. For t > i/eafc 
we observe a rapid decrease of P{t). The decay is still 
power-law but the exponent is larger than for t < U^ak 
though smaller than the classical result (7 « 1.5). This 
time scale can be utilized to study the interplay between 
quantum effects and classical anomalous diffusion. De- 
structive interference slows down the quantum motion 
however diffusion is not arrested due to the underlying 
classical anomalous (super)diffusion. 

Finally we address the similarities of our model with 
those of a generic Hamiltonian in a region of the phase 
space dominated by cantori. Qualitatively the effect of 
the trapping region is similar to that of cantori [1, [l3| • In 
both cases the particle remains in a small region of the 
phase space for a long time until eventually is released 
to the chaotic sea. In addition the dependence of the 
classical diffusion constant D ^ [11| is identical in 
both cases. In the quantum realm cantori [l^ [l^ [23| 
induces slow power-law decay in quantities such as the 
return probability [18,] or the eigenstates themselves [l^ . 
In addition cantori have been linked to fractional scaling 
with h [1^ and anomalous dynamical localization [2l| . 



All these features have also been observed in the model 
studied in this paper. These similarities strongly suggest 
that the 2 — 5 KR can be utilized as a toy model to study 
the effect of cantori in classical and quantum mechanics. 



IV. CONCLUSIONS 

We have studied the effect of classical anomalous diffu- 
sion on the quantum dynamics of atoms exposed to pairs 
of (5-kicks. Our result are generic, do not depend on the 
parameters Kt <^ 1 and r ^ 1 utilized. We have iden- 
tified a regime of quantum diffusion where the motion 
is slow down due to destructive interference but full dy- 
namical localization has not occurred yet. This has been 
related to the effect of the underlying classical trapping 
and releasing mechanism on the quantum dynamics. We 
have argued that the can be used as a simplified 

model to study the effect of cantori in classical and quan- 
tum mechanics. 
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